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We study the ordering kinetics in d = 2 ferromagnets which corresponds to populated neuron 
activities with long-ranged interactions, V (r) ~ r~ n associated with short-ranged interaction. We 
present the results from comprehensive Monte Carlo (MC) simulations for the nonconserved Ising 
model with n > 2. Our results of long-ranged neuron kinetics are consistent with the same dynamical 
behavior of short-ranged case (n > 4). The calculated characteristic length scale in long-ranged 
interaction is found to be n dependent (L(t) ~ t 1 /^" -2 )), whereas short-ranged interaction follows 
L(t) ~ t 1 / 2 law and approximately preserve universality in domain kinetics. Further, we did the 
comparative study of phase ordering near the critical temperature which follows different behaviours 
of domain ordering near and far critical temperature but follows universal scaling law. 


INTRODUCTION 


The spiking activity in complex neuron network in 
brain is dynamic (far from equilibrium) ]l|, exhibit 
nonequilibrium critical dynamics jH and the criticality 
in it has been used to characterize brain signals [dj. The 
reason could be when such system is quenched from a 
homogeneous phase to a broken-symmetry phase, it be¬ 
comes thermodynamically unstable. The subsequent far- 
from-equilibrium evolution of the system is characterized 
by the emergence and growth of domains enriched in 
the new equilibrium phases. This nonequilibrium evolu¬ 
tion, usually called kinetics of phase ordering or domain 
growth , has been the subject of much active investigation 
|4|. The domain morphology is quantified by the time 
dependence of the domain scale L(t), where t is the time 
after the quench. There is a good understanding of do¬ 
main growth kinetics in pure and isotropic systems with 
short-ranged interactions, where the domain scale shows 
a power-law behavior, L{t ) ~ t$ 0 [f>]. For the case with 
nonconserved order parameter, e.g., ordering of a ferro- 
magnet into up and down phases (spin-flip Glauber-Ising 
model 0 ), one has </> = 1/2 0 0 . On the other hand, 
for the case with conserved order parameter, e.g., phase 
separation of a binary ( AB ) mixture into A- and B-rich 
domains (spin-exchange Kawasaki-Ising model jlO|), we 
have 0=1/3 when growth is driven by diffusion 11, [12J ■ 
Apart from the domain growth laws, experimentalists 
are also interested in quantitative features of the domain 
morphologies. An important experimental quantity is the 
time-dependent correlation function C(r, t) or its Fourier 
transform, structure factor S(k, t) ( k being the wave vec¬ 
tor ii- Most of the studies are concentrated around 
the nearest-neighbor (short-range) inter-molecular inter¬ 
action. 

The neuron activities in brain at critical point are be¬ 
lieved to be effective for the long distance communication 
of the neurons [l3j because of coupling and variability to 


optimize information storage in the system 14] and dy¬ 
namic range of the system to response the input signal 
iTFj . The core of the paper focuses on long-range interac¬ 
tions which dynamically explains the rapid movement of 
the signal information inside the brain. The fast emer¬ 
gence of long-range interactions, mimic the rapid neu¬ 
ronal interactions in the brain. Further, we study the 
phase ordering dynamics in neurons with a specific inter¬ 
est to understand the role of the range of inter-neuron 
interactions. We address two important questions in this 
context via kinetic MC simulations: (a) What is the 
growth law for ordering phases of neurons? Is the growth 
law independent of the range of interaction? (b) What is 
the morphology for ordering phases of neurons, as mea¬ 
sured by the correlation function and structure factor? 
Is it comparable for all interaction range? We will be 
providing the answers to the above questions from our 
extensive MC simulations. 


NEURON ACTIVITY PATTERN MODEL 


Brain can be considered as a complex network of 
neurons which can be mapped onto ferromagnetic Ising 
model with two spin interactions 0 , where random fir¬ 
ing or non-firing of neuron can be represented by two 
states of a spin, s = +1 for firing and s = — 1 for rest or 
non-firing neurons 17]- Even though neuron activity pat¬ 
tern model is far from equilibrium dynamic model [18| . 
Ising model can be serve an excellent model to deal with 
critical phenomena of neuron activity pattern 0. The 
large number of local (short range) interaction of neurons 
tl9l ] and significant amount of global (long range) inter¬ 
action of neurons [0 are main basis of neurons commu¬ 
nication in brain network 01]. We consider the following 
Hamiltonian of two dimensional Ising system which in- 
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corporates long-ranged spin (neuron) model (LSM), 

H = - ^2 J(rij,n)siSj , Si = ±1, (1) 

<ij> 

where J is the coupling strength, n characterizes the 
range of the interaction, nj = |f) — fj\, and s* denotes 
the spin variable at site i. We consider two state spins: 
Si = +1 denotes an up-spin (active neurons) and Sj = — 1 
denotes a down-spin (passive neurons). We consider only 
a ferromagnetic case, where J > 0 always. The case 
where J can be both > 0 (ferromagnetic) and < 0 (an¬ 
tiferromagnetic) is relevant to spin glasses. We associate 
stochastic dynamics with the Ising model by placing it 
in contact with a heat bath. The appropriate dynam¬ 
ics for the phase ordering problem is spin-flip kinetics or 
Glauber kinetics. 



FIG. 1. Plot of < Si > vs. T for n= 2, 3, 4, 6, and 12 as 
indicated. The magnetization drops-off sharply near the crit¬ 
ical temperature ( T c ) and then vanishes to 0 in the disordered 
high-temperature phase. 


If we consider_ interacting neurons through slowly de¬ 
cay potentials 22], In order to capture thermodynamical 
parameters, one can define the following potential func¬ 
tion (obeying power law functional form) 


[23|, 24 


Uin) = lim — 
N—iOO N 


N 


J 


( 2 ) 

Here n = d + er = 2-|-cr [2f| for two dimensional system. 
For short-ranged interaction, a > 2 and the system size is 
not much important, whereas for long-ranged interaction, 
0 < a < 2 and it depends on the system size. Since, the 
size of the system is N, rescaling J —>• J/N to the Curie- 
Weiss model 24] and using Euler-McLaurin sum formula 
Fill for TV 1, equation © can be written as, 

Vn Vn 


U (n) = Inn J - 

iV-s-oo ' (~2 

X=1 y=l V*T 


1 


y 2 ) 


(n-2)/2 > 


rVN 

« lim J2 d(=2) / drg(r)r 3 ~ n , (3) 

N-> oo J l 

where, g(r) is the pair distribution function such that 
g(r) « 1 for r » 1. Then the equation © becomes, 

( 2 ln(TV) for n = 4, 

U(n ) « lim J i 7^4 (l - N 2 ~ n / 2 ) for n > 4, (4) 

N ^°° l 137774 TV 2- "/ 2 for 0 < n < 4 

From equation one can see that U(n) is finite for 
n > 4 when N —>• oo, and the asymptotic behaviour of 
finite critical temperature T c frit , 

TAn) » =» £ (-G) , (5) 

where, ks is Boltzmann constant. This shows that 
T c (n) oc 1/n for short-ranged potential (n > 4), whereas 
for long-ranged potential, T c depends on the size of the 
system N as well as n given by, 

(6) 

and T c diverges with system size. Similarly, one can also 
calculate other thermodynamical parameters such as in¬ 
ternal energy, entropy, free energy per particle (neuron) 
etc. at this asymptotic limit. 


DETAILS OF SIMULATION 


Since it is very difficult to obtain exact analytical so¬ 
lution of this problem, we straightforward implement a 
MC simulation of the Ising model with spin-flip kinetics 
to understand the behaviour. In a single step of MC dy¬ 
namics, we choose a spin at random in the lattice of dis¬ 
tribution of spins. The change in energy AH that would 
occur if the spin was flipped is computed with the step of 
acceptance or rejection based on Metropolis acceptance 
probability 27], |28| given by, 


f exp(— j3AH) if AH > 0, 
\ 1 if AH < 0. 


( 7 ) 


where, /? = (fc^T )” 1 denotes the inverse temperature. 
One Monte Carlo step (MCS) is completed when this 
algorithm is performed N times (where N is the total 
number of spins), regardless of whether the move is ac¬ 
cepted or rejected. All our simulations have been per¬ 
formed on a d = 2 lattice of size L 2 ( L s = 512) with peri¬ 
odic boundary conditions in both directions. The statis¬ 
tical quantities presented here (e.g., correlation function, 
structure factor) are obtained as averages over 10 inde¬ 
pendent runs. Each run starts with a randomly-mixed 
state with equal numbers of up (fired) and down (inac¬ 
tive) spins (neurons), which corresponds to a mean mag¬ 
netization m = (si) = 0 . 
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t = 100 t = 500 



FIG. 2. Evolution snapshots of domain coarsening for n = 2, 
3, and 6 quenched at T = 1 below T c at time t = 100, 500. The 
snapshots are obtained from a Monte Carlo (MC) simulation 
of ordering kinetics in ferromagnetic system. The details of 
the MC simulation are provided in the text. 


n=2 n=6 




FIG. 3. (a) Scaling plot of C(r,t) vs. r/L for a phase or¬ 
dering dynamics in d = 2 for n — 2. The data sets (for 
t = 100, 200, 500) collapse onto a single master curve, (b) 
Similar plot for n = 6. 


We study LSM for several values of n, namely 2, 3, 
4, 6 , and 12. The critical ordering temperatures, T c (n) 
for each n case is shown in Fig. [TJ where the charac¬ 
teristic behavior of spontaneous magnetization (< Si >) 
is plotted against temperature (T). As expected, T c (n) 
(dotted lines) increases with decreasing n as evident from 
equation ©• Above T c the spontaneous magnetization 


vanishes, whereas below T c it takes a nonzero value, in¬ 
ducing the typical behavior of a ferromagnet. Therefore, 
the physical properties of such systems and so its phase 
states depend on the value for the magnetization, the pa¬ 
rameter which is termed as order parameter : an ordered 
phase in which the spins are aligned appears when m/ 0 , 
while m = 0 implies a disordered (or symmetric) phase. 
Since, T c ’s for n > 4 are very close to each-other and 
hence, exhibit qualitatively similar behavior (explained 
shortly). We thus consider n < 4 cases for the long- 
ranged interaction. For each value of n, we cut-off the 
interaction at r c = (2.5) 6 / n to accelerate our simulation 
fH. We stress that the simulations are numerically very 
demanding for larger cut-offs. We compute several sta¬ 
tistical quantities to characterize the system. These are 
described as follows. 




FIG. 4. (a) Plot of C(r, t) vs. r/L at t = 500 for n = 2,3, 4,6. 
(b) Plot of S(k,t)L ~ 2 vs. kL, corresponding to the data sets 
in (a). The reasonably good data collapse shows that the 
scaling functions do not depend on the interaction range. The 
solid line denotes the OJK function in Eq. m & 


The domain coarsening is characterized by a growing 
time-dependent length scale L(t). The domain morphol¬ 
ogy does not change with time, apart from a scale factor. 
A direct consequence of the existence of a unique length 
scale is that the system exhibits a dynamical-scaling in 
the correlation function and structure factor. We com¬ 
pute the time-dependent correlation function: 


C (,Ti,rj;t) — {siSj) (si) (sj). (8) 

Here, the angular brackets denote an averaging over the 
independent initial ensemble and different noise realiza¬ 
tions. As the system is translationally invariant, the cor¬ 
relation function depends only onf= fj — f : 

C ( u, fj ■t) = C(fi,f i + f;t)=C (r, t). (9) 

Usually most experiments study the structure factor, 
which is the Fourier transform of the real-space corre¬ 
lation function: 

S(k,tj = j dfC{f,t)e itf! . (10) 

Since the system is isotropic, we can improve statistics 
by spherically averaging the correlation function and the 
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structure factor. The corresponding quantities are de¬ 
noted as C (r,t) and S(k,t), respectively. The correla¬ 
tion function and structure factor obey the dynamical 
scaling forms: 


C(r,t) =g[r/L{t )], 
S(k,t)=L(t) d f[kL(t)}. (11) 


Here, g(x) and f(p) are scaling functions; r is the sep¬ 
aration between two spatial points; k is the magnitude 
of the wave vector; and d is the system dimensionality. 
The characteristic domain size L(t ) is obtained as the dis¬ 
tance over which the correlation function decays to some 
fraction (say half) of its maximum value [C(r, t) = 1 
at r = 0]. There are several other suitable definitions 
for computing L(t), e.g., first zero-crossing of C(r,t ), in¬ 
verse of the first moment of S(k, t). In the scaling regime, 
all these definitions differ only by constant multiplicative 
factors (Sill. 



t 

FIG. 5. Time-dependence of the characteristic length scale 
L(t) for n = 2,3,4, 6, plotted on a log-log scale. The lines 
of slope 1/2 indicates the power-law growth regimes expected 
for phase ordering in d = 2 ferromagnetic system. 


NUMERICAL RESULTS 

In Fig. [2] we show the evolution snapshots obtained 
from our MC simulations for n = 2, 3, 6 with T = 1 
(< T c , see Fig. [TJ at t = 100, 500 MCS. At low tem¬ 
peratures, energetic effects are dominant and the system 
minimizes its energy by ordering the spins parallel to each 
other. In the absence of an external field (e.g., magnetic 
field, h = 0), the activated neuron (up-spin) and inacti¬ 
vated neuron (down-spin) states are equivalent. In the 
mean-field (MF) limit, i.e., n = 0, all the spins interact 


with each other and there is no spatial structure in the 
evolution morphology. For larger values of n, we see the 
emergence and growth of domains of up-spin (marked in 
black) and down-spin (unmarked). These domains inter¬ 
act and annihilate, resulting in coarsening of the char¬ 
acteristic length scale, and therefore, domain patterns 
at different times look statistically similar, apart from a 
global change of scale. The domain size at a fixed time 
(e.g., t = 500) is smaller for larger values of n. 

In Fig. [3l we show a scaling plot of the correlation func¬ 
tion, defined in Eq. ©. We plot C(r, t) as a function of 
the scaled distance r/L at three times, as indicated. Fig¬ 
ure ©a) corresponds to n = 2, and Fig. [31(b) shows data 
for n = 6. The data sets show a neat scaling collapse, 
confirming the validity of dynamical scaling. 


1=100 t = 500 



FIG. 6. Evolution snapshots of phase ordering systems in 
d = 2 for n = 2, 3, and 6. The system is quenched at T ~ T c . 
The MC simulations are described in the text. 

Let us next discuss whether the evolution morphology 
depends on the range of the interaction characterized by 
n. Fig. 2] shows a comparison of the scaling functions for 
four different n values (n = 2, 3, 4, 6) at a time t = 500, 
when the system is already in the scaling regime. In 
Fig. ©a), we plot the scaled correlation functions. The 
reasonably good data collapse suggests that the scaling 
functions do not depend on the interaction range. The 
solid line in Fig. ©a) denotes the analytical result due to 
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Ohta et al. (OJK) [33j, who studied ordering dynamics 
in a ferromagnet. The OJK function is 

C(r,t) = - sin- 1 (e- r2/i2 ). (12) 

7T 

(The corresponding result for the case with vector or¬ 
der parameter has been obtained by Bray and Puri [34|.) 
Our correlation-function data is in excellent agreement 
with the OJK function, showing that the phase ordering 
dynamics for n < 4 lie in the same dynamical univer¬ 
sality class as that for n > 4. In Fig. 0b), we plot the 
scaled structure factor [L~ 2 S(k,t) vs. kL] for the same 
time as in Fig. 0a). Again, the data sets collapse neatly 
onto a single master curve, confirming the scaling form 
in Eq. m- The scaling function is in excellent agree¬ 
ment with the corresponding OJK function. Notice that 
the structure factor, for large values of k , follows the 
well-known Porod’s law, S(k,t) ~ k~^ d+ ]k, which results 
from scattering off sharp interfaces [35|, [36|. The scaled 
correlation function and structure factor, in congruence 
with scale free behaviour of functional brain networks f371] 
depicts the universality of the interaction mechanism in 
both short and long range interactions in brain. 

In Fig.0 we turn our attention to the time-dependence 
of the domain size. We plot L(t) vs. t on a log-log scale 
for n = 2, 3, 4 and 6. Here, the data sets are consistent 
with the Cahn-Allen growth law, L(t) ~ f 1 / 2 — there is no 
sign of a crossover in the growth law at n = 4, as pre¬ 
dicted by Bray [38j. Bray has used the renormalization 


group (RG) approach to study ordering dynamics with 
long-ranged interactions of the form r d+a with 0 < <r < 2. 
In our case, d = 2 and a = n — 2. Bray argues that the 
long-ranged interactions are relevant for 0 < a < 2 or 
2 < n < 4, and irrelevant for n > 4. The corresponding 
growth law is 

f tVG-z) for 2 < n < 4, 

( f 1 / 2 for n > 4, 


m 


with possible logarithmic corrections. As we can see that 
our numerical results are not consistent with this pre¬ 
diction. The only difference as n is varied is that we 
have faster growth (higher prefactors) for smaller n, cor¬ 
responding to more long-ranged interactions. The fast 
dynamics of long-range interactions signifies the path of 
information processing and neuronal connections in the 
brain. The longer persistance of long-ranged neural con¬ 
nections could give sense to clustering behaviour of neural 
circuitry , specifically during learning of a specific task, 
new synaptic connections tend to form in vicinity of old 
connections related to that task [39;] making it more ro¬ 
bust. 

In Fig.0 we show the evolution of the order parameter 
(m) near critical temperature (T ~ T c ) from a disordered 
initial state (m = 0) for n = 2, 3, and 6 respectively. At 
higher temperature below T c we observe large fluctua¬ 
tions in the evolution patterns with very small global or¬ 
dering; instead of picking one of the up-spin, down-spin, 


or zero order parameter states, the system near T c is a 
kind of fractal blend of all three [ioj. However the cluster 
size is larger for smaller n. Recall that at T = 1 (<C T c ), 
thermal energy ( ksT ) of the system is low, thus spins 
try to obtain minimal energy by forming domains with a 
global ordering: m = +1 or -1. 




FIG. 7. (a) Plot of C(r,t) vs. r for n m 2, 3, 6 at t m 500 
MCS. (b) Sacling plot of C{r, t ) vs. r/L, corresponding to the 
data sets in (a). The solid line denotes the OJK function for 
the scaled correlation function in Eq. m- 

Finally, in Fig. 0a), we show the plot of correlation 
function [C(r, t) vs. r] corresponding to the evolution 
shown in Fig. 0 at t = 500. Note that the decay range 
of the correlation function is larger for smaller n. Fig¬ 
ure mb) shows the scaling plot of the data sets in (a). A 
reasonable data collapse confirms the dynamical scaling 
and clarifies that the system for each n belongs to the 
same dynamical universality class. 


SUMMARY AND DISCUSSION 

Let us conclude this paper with a summary and discus¬ 
sion of the results presented here. We study the effect of 
interaction range on the morphology of neuron activity 
pattern. The long-ranged and short-ranged interaction of 
neurons could be the main basis of how brain performs 
complicated functions at fundamental level. Neuron ac¬ 
tivities as well as wiring and rewiring of the neurons in 
the network subjected to heat bath depend on the range 
of interaction which are reflected in the dependence of 
critical temperature T c and magnetization on n. The do¬ 
main sizes of the neurons in short range interaction at 
far-lower critical temperature are smaller; some are iso¬ 
lated and numbers are more as compared to long ranged 
interaction for any time domain. However, the domain 
dynamics both in short and long ranged interaction sys¬ 
tem is quite different as compared to far critical temper¬ 
ature dynamics due to emergence of more randomness in 
the domain organization in the system. This leads to the 
change in domain growth laws of the neurons in short and 
long ranged interaction in the system. The correlation of 
neurons decays much faster in short ranged interaction as 
compared to long ranged, but it scales with r/L showing 
the universality of neuron interaction in brain. 
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Given the current focus on the biological network and 
their functionality, we hope that this paper will motivate 
fresh interest in the evolution dynamics and morphology 
of active and passive neurons. These kinetic processes 
play an important role in determining the functionality 
of brain. We emphasize that one can gain a good un¬ 
derstanding of the relevant neuron dynamics (wiring and 
rewiring inside the brain network) from simple coarse¬ 
grained models of the type discussed here. Convincing 
to the fact that re-learning help us to memorize things 
for longer duration. Thus it is an attempt to predict and 
an outlook to understand the functionaliies of the brain. 
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